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Abstract 

The pion scalar radius is given by (r|) = (B/ti") /^'^^2 dt 5s{t)/t'^ , with 

5s the phase of the scalar form factor. Below KK threshold, 5s ~ So, 5o being 
the isoscalar, S-wave tttt phase shift. Between KK threshold and t^^^ ~ 1.5 GeV 
I argued, in two previous letters, that one can approximate 5s ~ So, because 
inelasticity is small, compared with the errors. This gives (r|) = 0.75 ±0.07 fm^ 
and the value U = 5.4 ± 0.5 for the one-loop chiral perturbation theory constant, 
compared with the values given by Leutwyler and collaborators, (r|) = 0.61 ± 
0.04 fm^ and U = 4.4 ± 0.3. At high energy, t^^'^ > 1.5 GeV, I remarked that the 
value of Ss that follows from perturbative QCD agrees with my interpolation and 
disagrees with that of Leutwyler and collaborators. In a recent article, Caprini, 
Colangelo and Leutwyler claim that my estimate of the asymptotic phase 5s is 
incorrect as it neglects higher twist contributions. Here I remark that, when 
correctly calculated, higher twist contributions are likely negligible. I also show 
that chiral perturbation theory gives I4, — 6.60±0.43, compatible with my estimate 
but widely off the value U — 4.4 ± 0.3 of Leutwyler and collaborators. 
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1. Introduction. 



I will be concerned here mainly with four papers: they will be denoted by Yl,[^l Y2,[^l EYI"^! and CCL.W 
Fv{t), resp., Fs{t) will be the vector and scalar form factors of the pion; (r^) and (r|) the respective square 
radii. 5v{t) will be the phase of Fv{t) and Ss{t) the phase of Fs{t). Finally, Si{t), 6o{t) are the P and SO 
phase shifts for ww scattering. The pion scalar radius is given by 

{rl) = {6M / dt5s{t)/t'. 

For the vector radius one has a similar formula, 

/•oo 

{r^y) = (6/7r) / dt6v{t)/t\ 

J AMI 

Below KK threshold, 5s{t) = 5o{t) and 5v{t) = 5i{t). In fact, this second equahty probably holds as a good 
approximation up to t ~ 2 — 3 GeV^, where inelasticity begins to be important. 

Between KK threshold and t^/"^ ~ 1.5 GeV I argued, in Yl and Y2, that one can approximate 
~ (^0) because inelasticity is small compared with the errors. Using the experimental value for Sq, this 
gives (r|) = 0.75 ± 0.07 fm^ and the value f4 = 5.4 ±0.5 for the one-loop chiral perturbation theory (ch.p.t.) 
constant li, to be compared with the values given by Leutwyler and collaborators:'^'^! (r|) = 0.61 ± 0.04 fm^ 
and li = 4.4 ± 0.3. 

I also argued in Yl that at high energy, t > 2 GeV^ , the value of l)s that follows from perturbative 
QCD agrees with my interpolation and disagrees with that of Leutwyler and collaborators. In a recent 
article, Caprini, Colangclo and Leutwyler (CCL'^I) claim that my estimate of the asymptotic phase 5s is 
incorrect as it neglects higher twist contributions. Here I remark that, when correctly calculated, higher 
twist contributions are very likely negligible. I also show that ch.p.t. gives the number [4 = 6.6 ± 0.4, 
compatible within reasonable two loop corrections with my estimate; but widely off the value = AA ± 0.3 
of Leutwyler and collaborators. 

2. The asymptotic Fv{t), Fs{t) 



It is gcuerally assumed that one can approximate, at high energy, Fv{t) by its asymptotic form at leading 
twist (l.t.) the expression for which, neglecting quark and pion masses, has been known for a long time: 



t — *oo 



33 -2n/ -t\ogi-t/A^)' 



(For the definition of leading twist and twist three contributions to the pion form factor, see ref. 3). 

In 1983, Espriu and I calculated the twist three (t3) correction to this,'^! obtaining an infrared 
divergent result, 

^[log(_^/Mo)]^+^'^-; (2) 

t^oo (m„ + m.d) -t 

the m are the invariant quark masses, m(t) — ^mlog'^'" t/A^, and Mq is a cutoff mass (that may be taken 
equal to A or TOtt; it is irrelevant which, as the difference is subleading). We were not able to calculate the 
numerical constant C; for the other terms, cf. EY, especially the expressions for Bo{v) and the first formula 
in the second column of p. 189 in ref. 3. So, the full form factor would be 

the last because, at large t, the (t3) contribution dominates the (l.t.) one. Note that the (t3) contribution 
also dominates the (l.t.) one in the chiral limit, for all t. 
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When the paper was in the preprint stage, we received a letter from Brodsky and Lepage, pointing 
out that infrared divergent terms could not be treated as we had done; one should instead introduce a cutoff 
in ^ (the fraction of momentum of the quarks inside the pion) and treat the region near the endpoints as in 
refs. 6. If one does this, the twist three contribution becomes 



instead of (2), and the (l.t.) piece dominates: one recovers (1), for all of Fy- 

In EY^^] we were not fully convinced by the Brodsky-Lepage arguments [although we quoted them, 
including the result in Eq. (4) in the present paper], because a fit to data with (1) was much worse than a 
fit with (2) or (3), as we remarked in EY. Thus, in EY we concluded that "one really cannot calculate the 
pion form factor in perturbative QCD." 

However, some time later, Isgur and Llewellyn Smithl^l considered the partonic wave function of 
the pion (f>(^,t), for finite t. They concluded that the asymptotic limit is reached slowly and that, when 
reasonable (pi^jt) for finite t are used in the (l.t.) expression for Fy, the corresponding (l.t.) expression fits 
data well, even at subasymptotic energies. [An important remark is that the phase of Fy is not affected by 
the change in the wave function piece, which is real]. Because of this, in ref. 8 I used much milder words: 
I merely wrote that "There are, unfortunately, a number of snags" with the behavior (1), and proceeded 
to quote the EY, Brodsky-Lepage and Isgur and Llewellyn Smith!''! remarks, without drawing a definite 
conclusion. 

For Fs{t) the situation is very similar. To (l.t.) and leading order in perturbation theory, 

p(l■t■)/.^ ^ 5767r^ fUmu + rhdf 
S ^'t-oo33-2nf -i(log(-t/yl2)2dn. ' ^ ' 

as shown in Y2, Eq. (6.4). But, according to CCL, the (t3) contribution is^ 

^(t3).^. ^ 487rV2M2 log(-VMo) ^ 
^ t^oo 33 — 2n/ —t 

In this case, (6) would dominate over (5), both at large t and in the chiral limit, and thus, in particular at 
large t, we would have 

p m J. 576^2 f^{m^+md)^ ^Sir^f^M^ log(-t/Mo) 



iSn^f^M^ log(-t/Mo) 
*<x> 33 -2n/ -t 



(7) 



However, if one takes into account the Brodsky-Lepage remarks, the (t3) part should be calculated 
a la Brodsky Lepage Drell, would behave as t~^/'^, and hence would be negligible: one would recover 
the behaviour (5) given in Yl, for all of Fs- Note that the situation is identical for Fy and Fs; the 
difference between (l.t.) and (t3) contribution is due to the different pion wave functions involved, axial (p 
or pseudoscalar, (j)p, and these are the same for Fy and Fs- 

3. The asymptotic Sy{t), 5s{t) 



The phases of the form factors depend crucially on what happens to the (t3) contributions. If we believe the 
naive calculations in EY, or CCL, for the f i.e., that the behaviour of the form factors is given by (3) 

^ I am not sure of the correctness of the expression given in (6) for the (t3) piece; in EY we could not calculate 
analytically the twist-three wave function of the pion, <pp, hence we could not get the constant in Eq. (3), either. 
In fact we did not get constant asymptotic (pp but found it growing like a power of log t. I have the impression that 
CCL have forgotten that one has to diagonalizc the anomalous dimension matrix, i.e., that you have to calculate 
what we call Tnm in EY, which enters in the final result: see the detailed discussion in refs. 3, 8. However, since 
the conclusions are not altered by this, we will forget this problem here. 
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and (6), then, 

Svit)^^^n{l-ii + 2d„^)^-±^] [with(3)], 

(8) 



both phases approach the asymptotic value n from below. However, if we accept the Brodsky-Lepage 
arguments, one gets instead dominant (l.t.) and thus 

^>'<'>.-=„'{' + i^} [-"(1)1. 

fe(,),=,^.{l + 2,i„j-ljj-} [,ith(5)l, 

and the value tt is now approached from above. 

The last results, (9), were what I used in Yl, Y2. In fact, and because (as one can easily imagine) 
I was aware of the (t3) threat, I started by checking, in Yl, the behaviour (9) for Sv{t): here one has an 
independent evaluation of (r^), from the experimental pion form factor, 1^1 with which to compare. Looking 
at pp. 102, 103 in Yl, and correcting for an error there (which alters nothing numerically; see the Erratum 
in Yl), we find the following result: 

,2V f 0.434 fm^: with (9); ..q. 
1 0.303 fm^: with (8). ^ 

We have assumed the asymptotic region to begin at t = 3 GeV^ and taken 5v{t) — Si{t) below this. 
Experimentally, (rf.) = 0.432 ± 0.006 (the error here includes electromagnetic effects). The agreement 
of what one gets with (9) with experiment is very good; it only depends a little on the values of A or the 
matching point, that I took to be 1.7 GeV, provided we keep them reasonable. However, and this is more 
important, if using the naive (t3) result, given in Eq. (8), the value obtained is a disaster: below experiment 
by 30%. It is because of this that I took seriously, in Yl (and in Y2) the Brodsky-Lepage arguments and 
hence the (l.t.) evaluation of the phase of the form factors, Eq. (9). 

Nevertheless, I concede that the situation is not clear. While there exist a number of indications 
favouring the results of Yl, Y2 for (r|), one cannot definitely exclude the results of the evaluations of the 
Leutwyler group, as already remarked in Y2. One may therefore summarize the situation as follows: First of 
all, the experimental data for tttt scattering in the important region 2mK< t^'^ < 1.5 GeV are ambiguous. 
Different analyses produce mutually incompatible results, so the value for (rg) depends on which data one 
uses. If choosing the central value in the fit of Hyams et al.f^°l one gets (rg) = 0.61 ± 0.04 fm^ and hence 
= 4.4 ± 0.3 as in ref. 5. However, and as discussed in Y2, if one uses different sets of data (notably, taking 
into account information from tttt KK scattering) one finds a larger value, (rg ^) — 0.75 ± 0.07 fm^, 
leading to [4 = 5.4 ± 0.5. It could be possible, in principle, to discriminate between both solutions by using 
the asymptotic value for the phase of Fs{t) that follows from perturbative QCD; the result, however, is 
infrared divergent and thus, as happens for the twist three contribution to the electromagnetic form factor 
of the pion, a cut-off becomes necessary and the result is ambiguous. Therefore, one should probably accept 
a conservative estimate, covering both evaluations, 

r4 = 5±i. (11) 



4. li in chiral perturbation theory 



But this is not all. As I remarked in Yl, the number of low energy observables in tttt scattering, to one loop 
in ch.p.t., is such that one can determine the value of Z4 by fitting them, even if leaving Z4 free. 

In the present note I elaborate on this. To do so, I will use the values for scattering lengths and 
effective range parameters given in Table I. In units with the charged pion mass — 1, one has 



- 3 - 



-F. J. YNDURAIN- 





ABT 


ACGL [CGL] 


Exp. [PY] 


a,, 


U.ziy it U.UUo 


U.Z4U it U.UuU [U.zzU ± U.UUoJ 


U.zoU it U.UlU 




— U.U4Z ± U.UUi 


n HQf^ -I- n ni Q r n f\AAA -\- n nl 
— U.UoD ± U.Uio [— U.U444 ± U.UUiUJ 


— U.UOZ ± U.Uiz 




c\ o'7n _L r» ni 1 
U.z (9 ± U.Uii 


U.zoU ± U.Uz [U.z/o ± U.UUdJ 


n ofiQ _L n ni i 
U.ZDo ± U.Ull 


"o 


— U.U/D ^ u.uuz 


— U.U/y ^ U.UUD [ — U.UcSU ^ U.UUIJ 


n nsci -1- n ni 1 

— U.UoO ^ U.Uii 


10^ X ai 


37.8 ± 2.1 


37 ± 2 [37.9 ± 0.5] 


38.4 ± 0.8 


10^ X 61 


5.9 ± 1.2 


5.4 ±0.4 [5.67 ±0.13] 


4.75 ±0.16 


10^ X af^ 


22 ±4 


17ibl [17.5 ±0.3] 


18.70 ± 0.41 


10^ X af^ 


2.9 ± 1.0 


1.8 ±0.8 [1.70 ±0.13] 


2.78 ± 0.37 



CGL, ACGL: ref. 12, who apply dispersion relations together with experimental information at high energy and 
ch.p.t. (in CGL). ABT: ref. 13; these authors use models and chiral S'?7(3). Exp. [PY]: the experimental numbers as 
obtained in ref. 14 by fitting experimental data. The notation for the as and 6s is as in ref. 11. 

Table I 



These as and 6s can be described in terms of only four chiral coupling constants, Zj, i = 1, 2, 3, 4. 
In fact we will include in the fits the estimate of Gasser and Leutwylerl^^l for 1^ which gives 

Zs = 3 ± 2.5. (12) 

This has little influence the other parameters, but imposing (12) avoids unphysical values for Ts; otherwise 
the fits tend to give large negative values for this parameter, with very large errors. 

We turn to the determination of li from ch.p.t. We will use the one loop expressions for as and 6s 
in terms of U given in ref. 11, 

6(°)-^_/i + ^^r2/- +3r 131, ^',-1 

''^'^= 2887r3M./4{-^^+'2-^}- 
(for the D waves, see below) and we will consider several fitting strategies: 

i) Wc include in the fit the constraints (11, 12), as well as the experimental numbers (PY) for as and 6s 
given in Table I. Then, wc get a /d.o.f . — 8.8/(10 — 4) and the numbers 

h= -1.30 ±0.23, [2 = 6.04 ± 0.07, [3= 2.2 ± 2.5, [4 = 6.32 ± 0.61. (13a) 

The clear excess of the x^/d.o.f. over unity indicates that, at the level of accuracy of ref. 14, the two 
loop effects are not negligible. 
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ii) In fact, the errors in (13a) are purely nominal, as the fit depends on the way we have treated the higher 
order corrections. For example, if in the one loop expressions for a^^^ , b^^ , 



14407r3M^/4 



{[,+4r2-f }, 



_ ^ ff. r 1031. 

we replace the physical /n- by its value in the chiral limit, /, 



which is allowed since the difference is of higher order in ch.p.t., we find a /d.o.f. = 9.7/(10 — 4) and the 
central values 

h = -0.33, h = 4.56, [3 = 2.1, h = 6.93 : (13b) 
li and I2 are well outside the error bars given in (13a). 

iii) We can get safer error estimates, which take into account at least some of the uncertainty due to higher 
order effects, by averaging (13a) and (13b), weighted with the respective x^/d.o.f., and enlarging the error 
by including, as an extra error, the difference between the result of this operation and (13a). In this way we 
find an estimate that covers both (13a) and (13b): 

= -0.84 ±0.51, [2 = 5.34 ± 0.70, Zs = 2.2 ± 2.5, [4 = 6.61 ± 0.68. (13c) 



iv) The two loop corrections to nn scattering have also been calculated. [^"''1 It turns out that the scattering 
amplitude depends on six combinations of one and two loop coupling constants. Including two loop effects 
it is possible to give a virtually perfect fit to experimental tttt scattering, up to the F wave; but, unlike for 
the one loop case whore eight scattering lengths and effective range parameters arc fitted with only two Is. 
we now have six parameters for 14 or 16 observables (see below), so the agreement, although better, is less 
impressive than before. We find 

fi = -0.80 ±0.21, [2 = 5.5 ±0.35, [3 = 1-9 ±2.4, [4 = 6.60 ± 0.43. (13d) 

The numbers in (13d) have been obtained as follows: one fits the two loop expressions for scattering lengths 
and effective range parameters to the data deduced from experiment in ref. 14, for all waves up to the 
F-wave. One also includes in the fit the constraint given by Eqs. (13c) for the li. We then find a /d.o.f. = 
3.8/(14 - 10), if not including the F-wave in the fit, and a /d.o.f. = 7.9/(16 - 10) if including 4"^-' and 
63^'' in the fit. Averaging them and including their spread as an extra error (that can be attributed to three 
loop effects) we get (13d). 

This result in (13d), I4 = 6.60 ± 0.43, which is compatible with all others, is what we consider the 
more reliable one. 

v) The previous determinations include the constraint (11) in the fits. We could also fit leaving I4 completely 
free. In this case we find 

fi = -1.27 ±0.21, [2 = 6.04 ± 0.06, Z3 = 2.4±2.5, Z4 = 7.1±0.8, (14) 

with only the experimental errors in Table I, which is compatible with (13). We would find results similar 
to (13b, c, d) if including estimated higher order effects: the errors for I4 are slightly larger than before, and 
the central value is also only slightly larger. This confirms the result in (13d) for I4. 
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vi) Finally, we may also replace in the fits the "Experimental" PY data with those of ACGL or ABT in 
Table I. With only the errors there, and leaving I4 free, we get 

ri = -2.3±0.1, [3 = 6.20 ±0.05, [3 = 1.0 ±2.0, [4 = 6.6 ± 0.3, [CGL] 

/"i = -2.3 ±0.6, ^2 = 6.0 ±0.3, [3 = 2.9 ±2.5, [4 = 5.4 ±1.4, [ACGL] (15) 

fi = -1.32 ±0.55, fa =6.84 ±0.28, [3 = 1.1 ± 1.9, [4 = 6.71 ± 0.48. [ABT] 

The numbers that follow for I4, for all the fits, support the estimate in Yl, Y2,^ I4 — 5.4 ± 0.5, 
against that of Donoghue, Gasser and Leutwyler and Ananthanarayan et al.,'^' I4 = 4.4 ± 0.3. It is par- 
ticularly remarkable that the CGL evaluation of low energy parameters, which produces the parameters 
in the corresponding column in Table I, and which I used for the evaluation CGL in Eq. (15), was made 
imposing for I4 the value of the Leutwyler group, I4 = 4.4 ± 0.3: requiring consistency of the low energy 
parameters, however, pushes I4 to a value compatible with mine in Yl, Y2 and definitely incompatible with 
the I4 — 4.4 ± 0.3 input. In fact, and apart from the fit using the ACGLI^^J low energy parameters, whose 
large errors for I4 preclude drawing any clear conclusion from it, all other fits agree, in central value and 
error, with (13d). To get results for the value of I4 from the scalar radius of refs. 4, 5 consistent with what 
one gets from tttt scattering, one would require that the two loop corrections to the results of the Leutwyler 
group for (r|) be of 50% or more. 



5. Reggeistics 



Some comments are also made in ref. 4 about Reggeistics. I will only note here the surprising fact that, 
although a wealth of experimental data exist for high energy tttt cross sections, f'^^l which were included 
in the fits in refs. 17, Caprini, Colangelo and LeutwylerW do not use them to derive what they call "our 
estimate" (see Fig. 2 in ref. 4) and they do not compare what follows from said estimate with the tttt 
data. The interested reader may find more detailed discussions of the Reggeistics of the Leutwyler group in 
refs. 14, 17, 18. 
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